Spin Correlations in the Geometrically Frustrated Pyrochlore Tb2Mo2 07 
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We report neutron scattering studies of the spin correlations of the geometrically frustrated py- 
rochlore Tb2Mo207 using single crystal samples. This material undergoes a spin-freezing transition 
below Tg ~ 24 K, similar to Y2M02O7, and has little apparent chemical disorder. Diffuse elastic 
peaks are observed at low temperatures, indicating short-range ordering of the Tb moments in an 
arrangement where the Tb moments are slightly rotated from the preferred directions of the spin 
ice structure. In addition, a Q-independent signal is observed which likely originates from frozen, 
but completely uncorrelated, Tb moments. Inelastic measurements show the absence of sharp peaks 
due to crystal field excitations. These data show how the physics of the Tb sublattice responds to 
the glassy behavior of the Mo sublattice with the associated effects of lattice disorder. 
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Geometrically frustrated magnets can exhibit unusual 
phenomena at low temperatures, such as spin glass [l|, |2| , 
spin ice[3|, \4, \^, and spin liquid[6|, |7| states. The frus- 
trated pyrochlore compound Tb2Mo207 stands at an in- 
teresting crossroads, being a representative of both the 
Mo-based family i?2Mo207 {R = rare earth) and other 
Tb-based compounds Tb2Af207 (M = metal). As a 
function of the _R-site radius, the i?2Mo207 compounds 
exhibit a metal-insulator transition between ferromag- 
netic metal states to spin glass insulators.'®, 9i] The fer- 
romagnetic metals (such as Nd2Mo207 and Sm2Mo207) 
show an interesting interplay between the spin configura- 
tion and the transport properties, which can give rise to 
anomalous Hall and Ncrnst signals. lOlllII The spin glass 
insulators (such as Y2M02O7 and Tb2Mo207) have gen- 
erated much interest due to the apparent lack of chemical 
disorder and the relation to geometrical frustration. [l|, [3 
Outside of the Mo-family, pyrochlore compounds with 
Tb as the rare earth have been extensively studied be- 
cause the combination of exchange and dipolar interac- 
tions coupled with moderate spin anisotropy can lead to 
novel ground-states depending on the metal ion. For ex- 
ample, Tb2Ti207 has a spin liquid ground statej2]j while 
Tb2Sn207 exhibits magnetic order[l^, and both show 
persistent spin dynamics down to the lowest measured 
temperatures. 7', 13', 14] 

Tb2Mo2 07 crystallizes in the Fd3m cubic space group, 
and both the Tb and Mo sublattices form three di- 
mensional networks of corner sharing tetrahedra. Spin- 
freezing is observed below Tg ~ 24 K, and diffuse scatter- 
ing from the Tb moments has been observed with neutron 
scattering on powder samples. [llll5l.llq Remarkably, the 
spins keep fluctuating down to the lowest temperatures 
measured in /xSR studies [17|. Currently, the spin struc- 
ture of the short-range order is not known, nor has there 
been a precise determination of the correlation lengths. 
Also, not much is known about the possible effects of 
disorder, if present. In this paper, we present neutron 



scattering measurements on a single crystal sample which 
answer some of these questions and also uncover new fea- 
tures related to the interaction between the Tb and Mo 
sublattices. 

Single crystal samples of Tb2Mo2 07 were grown using 
the floating-zone technique in an image furnace. A mix- 
ture of Tb407 (99.97%) and M0O2 (99.99%) was thor- 
oughly ground and pressed into feed rods, and the crystal 
growth was carried out in an Ar atmosphere. The sam- 
ples were confirmed to be single phase Tb2Mo207 by x- 
ray and neutron diffraction. The magnetic susceptibility, 
shown in the inset of Fig. 1 (a) , is consistent with previous 
results on powders. The difference between the zero-field 
cooled (ZFC) and field cooled (FC) susceptibility indi- 
cate a spin- freezing temperature of Tg ~ 24 K. Fitting 
the high temperature susceptibility data (T > 100 K) 
to a Curie- Weiss law yields a Curie- Weiss temperature 
of Qcw — 12 K and an effective moment of ^ 9.5 /is. 
The size of the moment is consistent with large Tb'^^ 
moments [18| dominating the response. The nature of 
the magnetic couplings is not so clear; currently, the 
values of the Tb-Tb, Tb-Mo, and Mo-Mo interactions 
in Tb2Mo207 are not known. In Y2M02O7 (where Y 
is non-magnetic) , the Curie- Weiss temperature of about 
—200 K suggests that the Mo-Mo interaction is strongly 
antiferromagnetic.y Y2M02O7 has a spin-freezing tem- 
perature {Tg ~ 22.5 K) similar to that in Tb2Mo207; 
hence, it is reasonable to assume that the Mo sublattice 
is primarily responsible for the freezing transition. Neu- 
tron measurements on Y2M02O7 have shown that dy- 
namical short range order of the Mo sublattice develops 
below room temperature. |2J To better understand how 
the Tb''+ ions respond to the glassy behavior of the Mo 
sublattice, we have carried out a series of neutron scat- 
tering measurements. 

The neutron experiments were performed using a 0.49 
gram single crystal of Tb2Mo207 at the NIST Cen- 
ter for Neutron Research. The SPINS and BT9 triple 
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FIG. 1: (color online) Elastic scattering measured on SPINS 
spectrometer, (a) Intensity at the (002) diffuse magnetic peak 
position and the "background" signal measured at the (112) 
position as a function of temperature. Inset: Magnetic sus- 
ceptibility measured with H = 10 Oe along the [111] direction 
using a SQUID magnetometer. Scans through the (b) (002) 
position and the (c) (111) position at various temperatures. 
The horizontal bar indicates the instrumental resolution 

axis spectrometers were used with fixed final neutron 
energies of 4.5 meV and 14.7 meV, respectively. At 
SPINS, we employed the horizontally focused analyzer 
with collimator sequence 80 —sample— Be— radial— open. 
At BT9, a flat analyzer was used with coUimations 
40'-48'-PG-sample-40'-90'. The crystal was ori- 
ented in the {H H L)-zone and placed in a ^He cryostat. 

On SPINS, we observed the onset of diffuse elastic 
peaks upon cooling, and representative scans through the 
(002) and (111) positions are shown in Fig. 1(b) and (c). 
At T = 1.5 K, the peak widths are large, spanning a con- 
siderable fraction of the Brillioun zone, indicating short- 
range order. Since the magnitude of the Tb moment 
should be about 5 times larger than the Mo moment, we 
identify these peaks as arising from short-range order on 
the Tb sublattice. The intensity at the (002) diffuse peak 
position is plotted as a function of temperature in Fig. 
1(a), and the "background" signal measured at (112) is 
also plotted. The onset of the diffuse peak intensity oc- 
curs around T ^ 2Tg and increases smoothly upon cool- 
ing through Tg. This is similar to the temperature de- 
pendence of elastic scattering observed for Y2M02O7I3], 
indicating that the Tb correlations are coupled to the 
developing Mo correlations. 

In addition, we find a significant temperature- 
dependence of the Q-independent "background" signal 
in Tb2Mo207. Coohng from T = 109 K to around 50 K, 
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FIG. 2: (color online) (a) Correlation lengths of the short 
range order along different symmetry directions as a func- 
tion of temperature. Inset: Representative scans through the 
(111) position measured on the BT9 spectrometer, where the 
horizontal bar indicates the instrumental resolution, (b) Tem- 
perature dependence of the diffuse peak intensities. 

the intensity at (112) increases roughly linearly with de- 
creasing temperature. The Q-independence of the sig- 
nal can clearly be seen in the scans in Fig. 1(b) and (c). 
Typically, Q-independent background processes produce 
scattering which increases with increasing temperature, 
which is the opposite of what we observe. A likely ori- 
gin of this Q-independent scattering is the presence of 
frozen Tb moments which are completely uncorrelated. 
This identification is supported by the flattening of the 
temperature dependence below ~ 30 K, which is the 
same temperature at which the peaks due to correlated 
Tb moments begin to rise quickly. Such (5-independent 
scattering from Tb is unusual and may be caused by 
their coupling to the slowly fluctuating Mo moments 
as discussed further below. We flnd that both the dif- 
fuse peaks and the Q-independent signal are resolution- 
limited in energy. Although the energy resolution was 
narrow (Ai? ~ 0.12 meV), fluctuations at a rate below 
~ 0.02 THz would be indistinguishable from being static. 
Further elastic measurements were performed on BT9, 
and representative scans through the (111) peak are 
shown in the inset of Fig. 2(a). The solid lines indicate 
fits to Gaussian lineshapes plus a sloping background. 
The sharp nuclear Bragg peaks, denoted by the open 
symbols, were not included in the fits. Scans were per- 
formed along the il, H- and L-directions to determine the 
static correlation lengths as a function of temperature. 
Interpreting the Gaussian linewidths within a finite-size 
domain model, we calculated the correlation lengths ^ 
(or linear domain sizes) as plotted in Fig. 2(a). The 
correlations are not completely isotropic, with ^aiong [110] 
slightly larger than ^aiong [001] • Both correlation lengths 
increase modestly upon cooling and span a distance of 
roughly the size of a unit cell. Again, no dramatic 




FIG. 3: (color online) (a) Pattern of the elastic scattering 
intensity at T = 1.5 K measured on BT9 in the (HHL) zone. 
Data points due to scattering from nuclear Bragg peaks and 
Al powder lines (indicated by white regions) were removed 
in order to emphasize the diffuse magnetic scattering, (b) 
Calculated scattering pattern for short-range ordered regions 
with the spin- structure shown in the inset, which are domain- 
averaged as described in the text. A constant background of 
60 counts (which approximates the average background in the 
data) was added to the calculation. 



changes are observed at Tg. The fitted peak intensities 
are plotted in Fig. 2(b). The temperature dependence 
closely follows that of the SPINS data (also plotted) even 
though the energy resolution for the BT9 measurements 
is about 5 times broader. This further confirms that 
the diffuse peaks are associated with correlated moments 
which are static (fluctuating slower than ■~0.02 THz). 

An important issue is to understand the local struc- 
ture of the correlated Tb moments. In Fig. 3(a) we 
show the pattern of the elastic scattering intensity in the 
{HHL)-zone taken at T = 1.5 K on BT9. In previous 
powder measurements, two broad peaks were observed at 
IQI ~ 1 A^i and 2 A-i.[ll,[ll] Here, the single crystal 



data allow for a much more detailed examination of the 
intensities throughout reciprocal space. To first approx- 
imation, the peaks appear at positions consistent with 
a "(5=0 spin ice" arrangement (consisting of equivalent 
tetrahedra with moments along the local [111] anisotropy 
axes with a "two- in, two-out" arrangement). [J, |5| How- 
ever, a closer inspection of the intensities suggests devi- 
ations from the local spin ice structure. We calculated 
the magnetic scattering cross-section for various spin con- 
figurations on a small cluster of four Tb tetrahedra. A 
close match was achieved with an arrangement in which 
all tetrahedra in the cluster are identical and the Tb mo- 
ments are rotated away from the local [111] anisotropy 
axes of the spin-ice configuration by 14°. A depiction of 
this local order is shown in the inset of Fig. 3(b). The ro- 
tation of the Tb moments enhances the net ferromagnetic 
moment along the [001]-direction for each tetrahedra. 



Since there are six equivalent domains for this spin struc- 
ture, our calculation averages over all six magnetic do- 
mains. The uniform susceptibility data rule out a net fer- 
romagnetic moment at low temperatures, hence the fer- 
romagnetic moment from each cluster must cancel over- 
all, consistent with our domain-averaged cross-section. 
The calculated cross-section, plotted in Fig. 3(b), shows 
reasonably good agreement with the data. We estimate 
that the static moment associated with the diffuse peaks 
at r = 1.6 K is (Mrb) — 4.0(5) /is, significantly smaller 
than that expected for a free Tb^+ ion, which may partly 
be explained by crystal field effects. l^illSl This may also 
indicate that a substantial fraction of the Tb moment re- 
mains dynamic or is frozen without measurable spatial 
correlations. We comment more on the latter possibility 
below. Interestingly, our model for the local spin struc- 
ture is similar to the long-range order of Tb moments 
observed in diluted Tbi.8Lan.2Mo207 which is a ferro- 
magnet at low temperatures. fl6| 

Finally, we performed inelastic scattering measure- 
ments to investigate the Tb crystal field excitations, 
as well as possible collective excitations. In the struc- 
turally similar Tb2Ti207 and Tb2Sn207 compounds, a 
low-lying crystal field excitation has been observed near 
1.5 meV. This corresponds to an excitation from the 
ground state doublet to the first excited state. Surpris- 
ingly, in Tb2Mo207, we do not observe a sharp mode 
corresponding to this excitation. Figure 4(a) shows en- 
ergy scans taken at the (1.1,1.1,1) position on the SPINS 
spectrometer. Data taken at the (001) position (not 
shown) yield similar spectra. At T = 1.6 K, the data for 
positive energy transfers reveal a rather flat fluctuation 
spectrum without a clear peak in the range between 0.5 
meV and 4 meV. Warming to 20 K produces very little 
change; though further warming results in additional low- 
energy quasielastic scattering, consistent with previous 
observations [151] . The growth of this quasielastic scatter- 
ing appears to correlate with the diminishment of the Q- 
independent elastic signal. Inelastic scans were also per- 
formed at higher energy transfers (up to huj = 20 meV) 
on a 2.46 gram powder sample of Tb2Mo207 using the 
BT7 spectrometer, as shown in Fig. 4(b). No sharp peaks 
associated with crystal field excitations are observed; 
though, the low temperature spectrum shows broad fea- 
tures centered around 4 meV and 15 meV. In contrast, 
both Tb2Ti207 and Tb2Sn207 exhibit sharp crystal field 
excitations at about 10 meV and 15 meV.|7l. llQj 

The lack of sharp crystal field excitations for Tb'^+ 
in Tb2Mo207 is unexpected and requires better under- 
standing. One possible origin may be a static inhomoge- 
nous splitting of the non-Kramers Tb^+ doublet due 
to lattice strains. Such an effect has been observed in 
Pr2_2:Bia;Ru2 07 which has obvious chemical disorder. [20j 
While structural disorder has not yet been identified in 
Tb2Mo207, recent experiments on the Y2M02O7 com- 
pound reveal that there exists lattice disorder associ- 
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FIG. 4: (color online) Inelastic scattering data, (a) High- 
resolution measurements taken on SPINS at the (1.1, 1.1, 1) 
position in a single crystal, (b) Higher energy data taken on 
the BT7 spectrometer (with Ef — 14.7 meV and horizontally 
focused analyzer) on a powder sample at \Q\ = 3 A~^. Inset: 
Temperature dependence of the scattering at hco = 4.5 meV. 



ated with the Mo-Mo bond distances [2 ll |22| |. Since 
Tb2Mo207 has a similar spin-freezing transition, it is 
reasonable to expect similar levels of Mo sublattice disor- 
der. This would produce some degree of inhomogeneity 
in the Tb crystal field environments, thereby broadening 
the excitations. Another effect to consider is the inter- 
action of the slowly fluctuating Mo moments with the 
Tb moments. Freezing of the Mo moments manifestly 
breaks time-reversal symmetry, and therefore any cou- 
pling with the Tb sublattice would split the ground state 
Tb'^"'" doublet. The magnitude of the splitting would de- 
pend on the mean-field generated by the local configura- 
tion of Mo moments. Assuming the Tb^+ ground state 
doublet is primarily | ± 4) or | ± 5) as seen in related Tb 
pyrochloresflSl Il9| . then the mean-field at the Tb site 
due to the neighboring Mo would select a preferred local 
spin direction. As the Mo moments gradually become 
frozen in a disordered arrangement, the net effect on the 
Tb sublattice would be to select \ + Jz) or \ — Jz) in a 
spatially disordered way. Such a frozen and uncorrelated 
Tb configuration would appear as a Q-independent back- 
ground, consistent with our observations from the elastic 
scattering measurements. Indeed, the presence of static 
but uncorrelated Tb moments would naturally explain 
the relatively small value of {MTb) calculated from the 
diffuse peak intensities, since a significant fraction of the 
elastic signal would be Q-independent. 

In Fig. 4(b), the inelastic spectrum at T = 6 K of 



the powder sample exhibits enhanced scattering centered 
around 4 meV which quickly diminishes upon warm- 
ing to above ^^ 30 K (see the inset). The tempera- 
ture dependence closely follows that of the diffuse elastic 
peaks. This behavior suggests that the enhanced scatter- 
ing around 4 meV is derived from collective excitations of 
the Tb moments rather than single-ion physics. The en- 
ergy for this collective excitation should thus be related 
to a combination of the Tb-Tb and Tb-Mo interaction en- 
ergy scales. Further inelastic studies on the single crystal 
sample would be necessary to specify the magnitudes of 
the different interactions. 

In conclusion, neutron scattering studies on a single 
crystal of Tb2Mo207 reveal two components to the elas- 
tic scattering: a set of diffuse peaks plus a Q- independent 
signal. The short-range order of the Tb moments has 
been identified, and the observed correlation lengths are 
slightly anisotropic. The inelastic data indicate that the 
physics of the Tb sublattice is affected by disorder, most 
likely through lattice inhomogeneity and through cou- 
pling to the glassy Mo moments. The local clusters of 
correlated Tb spins that we observe may be the relevant 
unit in the slow dynamics found in this material 17[ and 
may connect to similar physics observed in other Tb- 
based pyrochlore systems as well 13|, llJ| . 
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